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Momentum distributions in time-dependent density functional theory: 
Product phase approximation for non-sequential double ionization 

in strong laser fields 
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We investigate the possibility to deduce momentum space properties from time-dependent density 
functional calculations. Electron and ion momentum distributions after double ionization of a 
model Helium atom in a strong few-cycle laser pulse are studied. We show that, in this case, the 
choice of suitable functionals for the observables is considerably more important than the choice 
of the correlation potential in the time-dependent Kohn-Sham equations. By comparison with 
the solution of the time-dependent Schrodinger equation, the insufficiency of functionals neglecting 
electron correlation is demonstrated. We construct a functional of the Kohn-Sham orbitals, which 
in principle yields the exact momentum distributions of the electrons and the ion. The product- 
phase approximation is introduced, which reduces the problem of approximating this functional 
significantly. 

PACS numbers: 31.15.Ew, 32.80.Rm 



I. INTRODUCTION 

Time-dependent density functional theory (TDDFT) 
[l| is a remarkably successful approach to the study 
of many-body systems in time-dependent external fields 
0. The essential statement of TDDFT is the same as 
that of the well-established ground state density func- 
tional theory (DFT) [|: all observables are, in princi- 
ple, functionals of the particle density alone. Since the 
latter is always a three-dimensional entity, independent 
of the number of particles involved, the computational 
cost of actual (TD)DFT calculations scales exponentially 
more favorable than the solution of the many-body (time- 
dependent) Schrodinger equation. 

In practice, almost all (TD)DFT calculations are per- 
formed using the (time-dependent) Kohn-Sham scheme 
[(TD)KS] (see, e.g., Q) where the density is calculated 
with the help of auxiliary, non-interacting particles mov- 
ing in an effective potential. The "art" of (TD)DFT is 
two-fold, namely finding sufficiently accurate approxima- 
tions to the density functionals of (i) the unknown effec- 
tive potential and (ii) the observables of interest. For- 
tunately, for many practical applications both items are 
uncritical Q ■ An example is the calculation of the opti- 
cal response of bio-molecules where even the simple lo- 
cal density approximation of the effective potential yields 
reasonable results, and the observable can be calculated 
from a known and explicit functional of the density (the 
time-dependent dipole). 

However, when it comes to the correlated motion of a 
few particles in a strongly driven system, TDDFT faces 
major challenges. In that respect, non-sequential dou- 
ble ionization (NSDI) serves as the "worst case" scenario 
for TDDFT. Theoretically, NSDI was addressed success- 
fully using the strong-field approximation (see^ e.g., [1] 
and references therein) and classical methods [a, @ . The 
widely accepted mechanism behind NSDI relies on the 



rescattering of the first electron with its parent ion, col- 
lisionally ionizing (or exciting) the second electron. 

In the recent publications Refs. 0, § significant 
progress was made in the treatment of NSDI within 
TDDFT as far as ionization yields are concerned. The 
latter display as a manifestation of the electron-electron 
correlation involved in NSDI the celebrated "knee" struc- 
ture in the double ionization yield, which was, until re- 
cently, not being reproduced within TDDFT. Reference 
[7| addressed issue (i) above (the effective potential) while 
Ref. [H focused on item (ii), the functional for the observ- 
able "double ionization". It was shown that (i) taking 
the derivative discontinuities at integer bound electron 
numbers into account and (ii) using an adiabatic approx- 
imation for the correlation function needed to calculate 
the double ionization probability, the NSDI "knee" can 
be reproduced. 

In our current work we turn to the much harder prob- 
lem of momentum distributions (or energy spectra 0]). 
In the NSDI regime the ion momentum spectra, as mea- 
sured in experiments employing "reaction microscopes" 
(see, e.g., 0, Ell), show a characteristic "double- hump" 
structure, i.e., maxima at non- vanishing ion momenta. 
The maxima at non-zero ion momenta are easy to under- 
stand within the rescattering scenario mentioned above: 
the first electron preferentially returns to the ion, colli- 
sionally ionizing the second electron, at times when the 
vector potential of the laser field is non-zero. Since the 
vector potential at the ionization time equals the final 
drift momentum at the detector, non-vanishing electron 
momenta (and, due to momentum conservation, non- 
vanishing ion momenta) are likely. In a TDKS treat- 
ment of NSDI in He starting from a spin-singlet state, 
the rescattering scenario is "hidden" in a single, spa- 
tial Kohn-Sham (KS) orbital. As we shall demonstrate, 
taking the auxiliary KS particles for real electrons and 
Fourier-transforming their position space product wave- 
function to momentum space leads to ion momentum 
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spectra in very poor agreement with the exact ones. A 
better approximation to calculate correlated electron mo- 
mentum spectra in the NSDI regime is required. With 
the present paper we aim at contributing to this goal by 
showing that item (ii) above, namely the construction 
of the functional for the observable, is the critical issue 
while (i) known effective potentials are sufficient, at least 
at the current level of accuracy. 

In Sec.|H] the model Helium system used to study the 
ionization process and the ensuing momentum distri- 
butions is introduced. In Sec. lIIII the method to cal- 
culate electron and ion momentum distributions is ex- 
plained. Results from the solution of the time-dependent 
Schrodinger equation (TDSE) in Sec. IIVI serve as a refer- 
ence for the results obtained using TDDFT in Sec.fVl 
The insufficiency of uncorrelated functionals to calcu- 
late electron and ion momentum distributions (Sec. IV A|) 
and the relative insignificance of the correlation poten- 
tial (Sec. IV B|) lead us to the construction of correlated 
functionals in Sec. IV CI In Sec. lVDl we introduce the 
product-phase approximation, which reduces the prob- 
lem of approximating the correlated functionals to that 
of approximating the exchange-correlation function. 

For consistency, we restrict ourselves to the presenta- 
tion of results for laser pulses with A = 780 nm and N = 3 
cycles. We stress, however, that the general conclusions 
drawn hold also for A = 614 nm, N = 3 and A = 780 nm, 
N = 4 laser pulses, as we have checked explicitly. 



II. MODEL SYSTEM 

A Helium atom exposed to linearly polarized laser 
pulses with N = 3 cycles and sin 2 -pulse envelopes is in- 
vestigated. The length of the pulses with a frequency 
of u> = 0.058 (corresponding to the experimentally used 
A = 780 nm) is T = 2Nit/uj, and the vector potential 
reads A(t) = A sin 2 t) sin(w t) for < t < T and zero 
otherwise (atomic units are used unless otherwise indi- 
cated). We use the dipole approximation, i.e., the spatial 
dependence of the laser field is neglected. The linear po- 
larization of the laser pulse thus allows us to describe the 
system by a one-dimensional model Helium atom with 
soft-core potentials for the Coulomb interactions. It is 
known that the essential features of the ionization pro- 
cess are described well by this model @, i, ED, El El El • 
Initially, the electrons are assumed to occupy the spin- 
singlet groundstate of Helium, and due to the neglect 
of magnetic effects in the dipole approximation the elec- 
trons stay in the spin-singlet state during the interaction 
with the laser pulse. Thus it is sufficient to study the 
spatial wavefunction, which has to be symmetric under 
exchange of the electrons. 

The TDSE idtip(xi,X2,t) — H(x\,X2,t) ip(xi,X2,t) is 
solved for laser pulses with different effective peak inten- 
sities / = 1(A). A trivial gauge-transformation cancels 
the purely time-dependent A 2 -term and yields the Hamil- 



tonian 

S= £ {-\d 2 Xi + V(x i ,t))+W(\x 1 -x 2 \), (I) 

8=1,2 

with H = H(xi,X2,t). The external potential is 
V(x, t) = —iA(t) d x — 2/Va; 2 + e C n, the electron-electron 
interaction potential is given by W(x) — 1/ ^Jx 2 + e GC . 
The soft-core parameters e cn and e ee are chosen to yield 
the correct ionization potentials. Reproducing the ion- 
ization potential of He + , = 2.0 in a corresponding 
model He + ion fixes e cn — 0.5. The choice e cc = 0.329 
yields the ionization potential of Helium, — 0.904. 
All results presented in this work are qualitatively insen- 
sitive to the precise values of the soft-core parameters. 

As the two electrons constitute a spin-singlet state 
for all times they are described by the same KS or- 
bital. Therefore, in a TDDFT treatment, we have 
only one time-dependent Kohn-Sham equation (TDKSE) 
id t 4>(x,t) = H KS (x,t) 4>(x,t) with the Hamiltonian 

H KS (xA) = ~d 2 x + V(x,t)+v hx (x,t)+v c (x,t). (2) 

The Hartree-exchange potential Uhx — Vh + v x follows as 
Vbx (x,t) = \ $ dx'n(x',t)/yj(x -x') 2 + iff. We have 
used the exact exchange term for Helium v x (x,t) — 
— v\ i (x,t)/2, which is local as both electrons are de- 
scribed by the same orbital. 

Setting v c — yields, in the special case of the 
Helium atom or He-like ions, an identical description 
as the time-dependent Hartree-Fock (TDHF) treatment 
(due to the locality of w x ). The LK05 potential v^ K05 
[3] takes into account the discontinuous change in the 
correlation potential when the number of bound elec- 
trons N(t) = J_ a dxn(x,t) passes integer numbers, 
v?™(x,t) = [B(t)/(l + exp[C(B(t)-2)])-l]vte(x,t), 
where C is a sufficiently large constant (we set C = 50) 
and B(t) = No/N(t). In order to encompass all bound 
states the parameter a is chosen as a — 6 a.u. throughout 
this work, results being insensitive to the precise value of 
a. We use e^f = 0.343 in the Hartree-exchange potential 
i>h x to acquire I^ r> = 0.904 for the model Helium atom. 
The TDSE and TDKSE are solved by a split-operator 
time propagator on a numerical grid (see, e.g., [15[ and 
references therein). 

Along the lines of Ref. Q we construct from the 
TDSE solution an exact KS orbital (EKSO). The 
Schrodinger solution gives the exact density of our 
model Helium atom n(x,t) — 2 j dx2 \ip(x, x 2 , t)\ 2 = 
2 J dx\ \ip(%i, x, t)\ 2 and the exact probability current 
j(x,t). From the equality of the exact and KS 
currents in the case of a one-dimensional system, 
the phase of the EKSO is determined as i?(a;,t) = 
J oo dx'j(x',t)/n(x' 1 t) + a(t). The unknown purely 
time-dependent phase factor a(t) does not affect the re- 
sults presented in this work and is therefore set to zero. 
The EKSO cj)(x,t) = y/n(x, t)/ 2e l ^ x ^ is thus identical 
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to the orbital a TDDFT calculation with the exact cor- 
relation potential v c would yield via the TDKS scheme. 
The EKSO allows us to separate the challenges facing 
TDDFT calculations (cf. Sec.|T|: finding (i) a suitable 
approximation of v c (where it serves as a reference for 
the resulting orbital) and (ii) appropriate functionals for 
observables (where it is the exact input). 

III. MOMENTUM DENSITIES 

We partition the two-electron space and associate 
with single ionization the area A(He + ) = {(xi,X2) | 
\xi\ > a, \xj^i\ < a V i,j G {1,2}} and with dou- 
ble ionization the area _4(Hc 2+ ) = {(xi, x%) \ \x\\ > 
a,\x2\ > a}. Integrating \ij){x\, X2, t)\ 2 over these ar- 
eas then yields the respective ionization probabilities, 
with the double ionization probability given by P 2+ (t) = 
J/^4( Hc 2+) dxi dx2 \^(xi,X2,t)\ 2 . This scheme to deter- 
mine ionization probabilities from the two-electron wave- 
function has been successfully used in numerous similar 
calculations [H H [3 . 

The wavefunction ip(xi,X2,t) can be de- 
scribed equivalently in momentum space by 
its Fourier transform (27r) ip(k\, k%,t) = 
J dxi f dx 2 ^(xi,X2,t)e- i(klXl+k2X2 \ As the wave- 
function in momentum space is normalized to one, 
the pair density in momentum space is given by 
p{k!,k 2 ,t) = 2\^{k 1 ,k 2 ,t)\ 2 . 

At times < t < T during the laser pulse the velocity 
of the electrons is actually given by ±i(t) = ki(t) + A(t), 
i.e., the sum of the canonical momentum fcj and the value 
of the vector potential at the respective time. In this 
work we investigate properties of the system at t = T 
after the laser pulse. As A(T) — 0, canonical momenta k 
and drift momenta are identical. 

We are interested mainly in the double ionization pro- 
cess and thus Fourier transform only the wavefunction 
in the area >4(He 2+ ) associated with double ionization. 
The resulting sharp step at the boundary of A (He 2+ ) at 
\xi\ = a, \xj^i\ > a, with i,j € {1,2}, is a potential 
source of artifacts when Fourier transformed. Hence, a 
smoothing function f(xi,X2) — Yii=i + e~° \ x i~ a \ 

is introduced. The factor c has to be of the order of 
one, in this work we choose c = 1.25. The smoothing 
function is constructed so that J J dxidx 2 f 2 (x\,X2)b = 
IIa (Hc 2+ ) dxi dx2 b for a constant b. This condi- 
tion ensures that the wavefunction ip( 2+ \xi, X2, t) — 
f(xi,X2)ip(xi,X2,t) gives to a good approximation the 
same double ionization probability as the original wave- 
function, i.e., that J dxi J dx2 f 2 (xi, X2)\tp{xi, x 2 , t)\ ~ 
P 2+ . The correlated wavefunction of the electrons freed 
in double ionization in momentum space is thus calcu- 
lated as 

(27T) V (2+) (fcl,fc2,t) = 

J d Xl J dx2^ {2+) (x 1 ,X2,t)e-' l{klxl+k2X2) . (3) 



This approach is equivalent to projecting out the states 
corresponding to single and no ionization and is known 
to lead to accurate momentum distributions [T3 |. 

From the wavefunction we construct the momentum 
pair density of the electrons freed in double ionization 

p ( - 2+ \k l ,k 2 ,t) =2\^ 2+ \k ll k 2l t)\ 2 . (4) 

The probability to find at time t an electron freed in dou- 
ble ionization with momentum k% in dfci and an electron 
with fc 2 in dk 2 is then p( 2+ \ki, k 2 ,t) dfcid/c2- 

In experiments, it is easier to measure the momentum 
of the He 2+ ion fc Ion after double ionization instead of in- 
dividual electron momenta. As the total photon momen- 
tum involved is negligibly small, this provides informa- 
tion about the sum of the electron momenta via momen- 
tum conservation k\ + k 2 = — k lon . The ion momentum 
density then follows from the momentum pair density of 
the electrons freed in double ionization (j4|) as 

n[ 2 J ] (k Ion ,t) = I J dkp^ 2+ \-k lon -k,k,t) 

= ~ J dkp( 2+ \k,-k lon -k,t), (5) 

due to the symmetry of the electron momentum pair 
density. The factor 1/2 ensures the correct normaliza- 
tion since the system consists of only one ion but two 
electrons. The ion momentum density n[ 2 ^(fc Ion , t) dk larL 
gives the probability to find at time t the He 2+ ion with 
momentum fc Ion in dfc Ion . 

IV. MOMENTUM DISTRIBUTIONS FROM 
THE TDSE 

From the numerical solution of the TDSE we obtain 
ip(x\, X2, T) after the interaction with the laser pulse. In 
the left hand side of Fig.[T] the momentum pair density 
of the electrons freed in double ionization, as calculated 
from Eq. ([J), is shown. 

For all but the highest intensity depicted, electrons 
have the highest probability to move at different veloc- 
ities I fci| ^ \k 2 \ (±i(T) = ki(T) since A(T) — 0, cf. dis- 
cussion in Sec. IIIip but in the same direction (sgn(fci) = 
sgn(fc2)). Depending on the laser intensity the proba- 
bility for the double ionization process is highest at dif- 
ferent half-cycles of the laser pulse, i.e., different signs 
of the vector potential. Therefore, the favored direction 
in which the electrons leave the atom varies with inten- 
sity. NSDI can be understood by a recollision mechanism 
where one electron returns to the He + ion and frees the 
second electron (see, e.g., jH). The results of the TDSE 
then imply that both electrons leave the atom in the same 
direction but due to Coulomb repulsion their velocities 
differ, in accordance with earlier results for a longer laser 
pulse [HI]. 

The "butterfly" shape of the momentum pair density of 
the electrons freed in double ionization as shown in Fig.[T] 
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is evidence that it is highly correlated, as it cannot be 
reproduced by multiplying two orbitals for the respective 
electrons. 

For I = 6.96 x 10 15 W/cm 2 both electrons have the 
highest probability to leave the atom in the same direc- 
tion with similar velocities k\ s» k% . This can only be the 
case when the Coulomb repulsion between the electrons 
is weak, i.e., when they are removed sequentially, result- 
ing in a large spatial separation. The final non- vanishing 
velocities are due to the high intensity of the laser pulse, 
which ionizes the atom so rapidly that A(t) ^ when the 
first electron is freed. The grid-like structure typical for 
a product wavefunction is seen, the electron correlation 
being weak. 

From the momentum pair density of the electrons freed 
in double ionization p^ 2+ '(ki,k2,T) (j4|) we calculate the 
ion momentum density (k\, k2,T) ([5|). For different 
effective peak intensities the density of the ion momen- 
tum is depicted in Fig. [2] It exhibits peaks at non-zero 
momenta. As explained in Sec. [J these are typical for rec- 
ollision processes when the first freed electron recollides 
close to the maximum of the vector potential, i.e., when 
\A(t)\ w A. Hence, the sum of the momenta of both elec- 
trons is non-zero, and, by momentum-conservation, this 
holds for the ion momentum as well [161 ]. 

For an infinitely long laser pulse of laser period T/N, 
H(t + T/N) = H(t) holds while this symmetry is bro- 
ken in the case of few-cycle laser pulses. Hence, with 
respect to the dislodged electrons there is no spatial in- 
version symmetry, leading to asymmetric ion momentum 
distributions 17 , HH, Ell ■ This effect is clearly seen in 
Fig. [2 For the three lowest intensities a process with 
fc Ion > dominates while with increasing intensities pro- 
cesses with fc Ion < become more likely. In addition, a 
central peak gets more and more pronounced, showing 
that the relative probability of sequential double ioniza- 
tion increases. The fact that the peak is not centered 
around fc Ion = for 7 = 6.96 x 10 15 W/cm 2 is again due 
to the high intensity and the short duration of the laser 
pulse, as explained above. 



V. MOMENTUM DISTRIBUTIONS FROM 
TDDFT 

DFT can be formulated in momentum space (see, e.g., 
(20|). and this seems to be the obvious path to follow 
when one is interested in the calculation of momentum 
spectra. However, momentum space DFT lacks the "uni- 
versality" feature of the Hohenberg-Kohn theorem 
meaning that each system under study requires a dif- 
ferent momentum space effective potential — an entirely 
unattractive feature. We therefore prefer to make the 
"detour" via standard, universal, position space TDDFT. 
In the case of single ionization, a straightforward cal- 
culation of the momentum or energy spectrum from 
the Fourier-transformed valence KS orbital may be a 
good approximation (see, e.g., the approach followed in 



Ref. [ll|). Instead, it is less obvious how to determine 
correlated momentum spectra from position space TDKS 
orbitals. 

As explained in the Introduction, determining momen- 
tum pair densities and ion momentum densities from a 
TDDFT approach faces two challenges: The first is to 
find an approximate correlation-potential v c in the TD- 
KSE to reproduce the exact density n(x, t) with sufficient 
accuracy. The second, more difficult one, amounts to as- 
sign a suitable functional of the density to the respective 
observable. As both the ion momentum density and the 
momentum pair density (via their probability interpreta- 
tions, cf. Sec. lIII[) are observables, the Runge-Gross the- 
orem assures that functionals of the density alone exist 

Q. 



A. Uncorrelated functionals 

Treating the KS orbital as if it were a one- 
electron wavefunction yields a product wavefunction 
4>{x\, t) (j>(x%, t). This is the same assumption frequently 
made to derive uncorrelated ionization probability func- 
tionals (see Ref. [|| and references therein). 

The Fourier transformed KS orbital for \x\ > a, i.e., 
with the bound states projected out (see Sec. lIIip is 



2 7T0 (+) (fc,f) = / dxf(x)(p(x,t)e- 



(6) 



with f{x) = l/\/l + c~ c \ x ~ a \ the one-dimensional 
smoothing function equivalent to the smoothing func- 
tion used in Sec. lIIII Calculating the momentum pair 
density ^ and the ion momentum density ([5]) from the 
product wavefunction gives the uncorrelated functional 
for the momentum pair density of the electrons freed in 
double ionization 

p< 2+ ) (ki,k 2 ,t) = 2\ <t> {+) (to , t) <M +) (h , t) | 2 (7) 

and the uncorrelated functional for the ion momentum 
density of He 2+ 

n[^\k lm ,t) = /dfc|^ + >(-fc Ion -MH (+) (M)| 2 - (8) 



Equations ([7]) and ([8]) are not functionals of the den- 
sity alone but due to the Fourier transformation they are 
dependent on the density and on the phase of the KS 
orbital. 

The momentum pair density at t — T, as calculated 
from the uncorrelated functional ||7J) using the EKSO, 
is depicted in the right part of Fig.Q] for A = 780 nm, 
N = 3-cycle laser pulses with different intensities. Com- 
parison with the left hand side showing the momentum 
pair density calculated from the correlated Schrodinger 
wavefunction ip{x\ ,X2,T) confirms that only for the high- 
est intensity a product wavefunction approach is reason- 
able. For lower intensities the uncorrelated functional 
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FIG. 1: Contour plots of the momentum pair density 
p 2+ (ki, lt2, T) of the electrons freed in double ionization. Re- 
sults calculated from the uncorrelated functional (|7j) using 
the EKSO (right hand side) are compared to the TDSE (left 
hand side) solution. Momentum pair densities for A = 780nm, 
N — 3-cycle laser pulses with different effective peak intensi- 
ties are shown. 
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FIG. 2: Ion momentum density of the model He ion after 
interaction with A = 780 nm, N = 3-cycle laser pulses with dif- 
ferent effective peak intensities. The density calculated using 
the EKSO in the uncorrelated functional © is compared to 
results from the TDSE. 



B. The role of the correlation potential 



for the momentum pair density does not exhibit the 
typical "butterfly" -shaped correlation structures of the 
Schrodinger solution. Instead, the grid-like structure typ- 
ical for a product wavefunction is clearly visible. 

For the same system we calculate from Eq. (jHJ the ion 
momentum density using the EKSO. In Fig. [2] the He 2+ 
ion momentum density is compared to the results from 
the TDSE, which are scaled to enable the comparison of 
qualitative features. The different values of the integrals 
over the ion momentum densities are due to the differ- 
ent double ionization probabilities, as can be seen from 
J dk lon n^\k Ion ,t) ~ P 2+ , which follows from Eq. © 
(see Ref. Q and references therein for a discussion of this 
particular problem). Apart from the highest intensity 
the density is centered around a central peak at k lan w 0. 
This is evidence that correlations, which are not included 
in the uncorrelated functionals for the observables, are 
responsible for the distinct peaks of the ion momentum 
density at non-zero momenta. This result is consistent 
with the analysis of the results of the TDSE (Sec. lIV|) . 
which attributes the peaks at k Ion ^ to electron rescat- 
tering, i.e., to an interaction between the electrons. For 
the highest intensity shown in Fig.[2j sequential double 
ionization becomes dominant (cf. Sec. lIV[) . so that the de- 
scription using the EKSO in the uncorrelated functional 
reproduces the ion momentum density reasonably well. 



To underline the importance of the functional for the 
ion momentum density we use the correlation potentials 
v c = (TDHF) and i£ KQ5 (LK05) in the TDKSE for our 
model He atom interacting with the A = 780 nm, N = 3- 
cycle laser pulses (cf. Sec.lH|). 

In Fig.[3]the ion momentum densities obtained from us- 
ing the respective orbitals in the uncorrelated functional 
for the ion momentum density 1(5)) are compared to the 
results with the EKSO, i.e., the orbital which the exact v c 
would yield. For the TDHF approach, results are similar 
to the results using the LK05-potential. Both approx- 
imations lead to uncorrelated ion momentum densities 
which are close in qualitative terms to the EKSO results. 
Only at the highest intensity I = 6.96 x 10 15 W/cm 2 they 
exhibit a single peak at k Ion > and not, as the EKSO 
solution, at k lon < 0. In this intensity regime purely 
sequential double ionization dominates, pointing to pos- 
sible shortcomings in the description of this process with 
both correlation potentials. 

As the general deficiencies of the uncorrelated func- 
tional described in the previous paragraph are entirely 
due to the functional for the observable, these results 
demonstrate the relative unimportance of the choice of 
the correlation potential in the TDKSE for the observ- 
ables of interest in this work. 
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FIG. 3: Ion momentum density of the model He ion after 
interaction with A = 780 nm, N = 3-cycle laser pulses with dif- 
ferent effective peak intensities. The densities are calculated 
from the uncorrelated functional © using the EKSO and the 
orbitals obtained with v c = (TDHF) and ^ K05 (LK05). 



C. Towards correlated functionals 

In polar representation, the solution of the TDSE is 
written as ip(x 1 ,x 2 ,t) — ^ p(x\,x 2 ,t) / 2 e ^<f(xi,x2,t) anc j 

the KS orbital as (j>(x,t) = v / n(x,t)/2e ii>ix ' t \ We define 
a time-dependent complex exchange-correlation function 



n(xi,x 2 ,t) 



tp(x 1 ,x 2 ,t) 



V2<f>(xi,t)<f>(x 2 ,t) 

/g^ e i[v>(x 1 ,X2,t)-tf(x 1 ,t)--d(x 2 ,t)] ^ 



with the time-dependent exchange-correlation func- 
tion g xc = g xc (xi,x 2 ,t) given by g xc (x 1 ,x 2 ,t) = 
p(xi, x 2 , t)/ n(xi, t) n(x 2 , t). Approximations to g xc = 
\k\ 2 have been used to construct correlated ionization 
probability functionals [1, [22J. Note that while g xc is 
an observable (and thus a functional of only the density 
exists), the complex- valued k is not an observable. Using 
Eq. ^ to express the correlated wavefunction ip(xi, x 2 ,t) 
in terms of the KS orbitals and the complex exchange- 
correlation function, Eq. (H|) gives the correlated func- 
tional for the momentum pair density of the electrons 
freed in double ionization 



dxi / dx 2 k(xi, x 2 , t) 



x ^ + \x 1 ,t)^ + \x 2 ,t)e-^ klXl+k2X ^ (10) 

with <p + > {x ,t) — f(x) 4>{x, t). The correlated ion mo- 
mentum density is calculated by using the correlated mo- 



mentum pair density in Eq. ([5]). We thus have exact mo- 
mentum distribution functionals, which depend only on 
the complex exchange-correlation function k and the KS 
orbital </>. 

The complex exchange-correlation function n in turn 
depends on the pair density and the phase of the 
Schrodinger solution ip(xi,x 2 ,t). In order to derive mo- 
mentum space properties for more complex atoms than 
Helium from the KS orbitals directly through expressions 
like Eq. (|10p . it is inevitable to approximate k. However, 
this is challenging since, due to the Fourier-integrals in 
Eq. pTJj) . the complex exchange-correlation function has 
to be approximated in all „4(He 2+ ) (and not just for the 
bound electrons, as in the calculation of ionization prob- 
abilities nm). 



D. Product phase approximation 

The necessary approximation of the complex exchange- 
correlation function k ([9]) consists of approximating 
g-Kc{x\,x 2 ,t) and the phase-difference ip{x\,x 2 ,t) — 

^{X!,t) - §{x2,t). 

Addressing the second part, the easiest approximation 
follows from the assumption that the difference of the 
sum of the phases of the KS orbitals and the phase of 
the correlated wavefunction can be neglected when cal- 
culating momentum distributions, i.e., we set 



ip(xi,X 2 ,t) = r?(xi,t) +-&(x 2 ,t). 



(11) 



Since $(a;,i) is the phase of the KS orbital we denote 
this approach as the product phase (PP) approximation, 
which yields 



K PF (xi,x 2 ,t) = y/g xc (xi,x 2 ,t) . 



(12) 



It is noteworthy that knowledge of the exact k pp thus suf- 
fices to calculate the exact double ionization probabilities 
from the EKSO. 

We calculate the ion momentum density using Eq. (fT2|) 
in Eq. dTUJ) and in Eq. ©. Employing the EKSO, the 
ion momentum densities shown in Fig. 3] for A = 780 nm, 
N = 3-cycle laser pulses with different intensities are 
obtained. The results from the TDSE are depicted as 
well. For comparison of the qualitative features, they 
are scaled, although the integrals over the ion momen- 
tum densities are equal in both cases (note that the PP 
approximation returns the exact double ionization proba- 
bilities). A generally good qualitative agreement with the 
Schrodinger solution is acquired. The asymmetric struc- 
ture and distinct peaks are reproduced. For intensities 
where NSDI is strongest, the quantitative agreement is 
least convincing. Although the PP approximation does 
not reproduce the exact fc Ion positions of the peaks, it 
modifies the uncorrelated functionals in a way which al- 
lows to deduce information about the underlying double 
ionization processes at the different intensities. We can 
therefore conclude that the difference between the phase 
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FIG. 4: Ion momentum density of the model He 2+ ion calcu- 
lated from the correlated functionals in the PP approximation 
using the EKSO. Results for A = 780 nm, N = 3-cycle laser 
pulses with different effective peak intensities are compared 
to the ion momentum density obtained from the TDSE. 
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FIG. 5: Contour plots of the exchange-correlation function 
gxc(xi,X2,t) for two effective peak intensities of A = 780nm, 
N = 3 cycle laser pulses as acquired from the solution of the 
TDSE. For clarity values larger than 10 are shown as 10. 



of the correlated wavefunction and a product wavefunc- 
tion is not as important for reproducing the structure 
of the ion momentum density as is the correlation given 
by g xc (xi, X2,t). This conclusion was verified by setting 
g xc = 1 in Eq. §§§ and using the exact phases in Eq. (fTU)) , 
which did not yield the peaks present in the Schrodinger 
solution. Using LK05 orbitals in the PP approximation 
also reproduces distinct peaks while the general agree- 
ment with the Schrodinger ion momentum density is not 
as good as for the EKSOs. 

The contour plots of the momentum pair density of the 



electrons freed in double ionization p^ 2+ \k\,k2,t) calcu- 
lated from the correlated functional in the PP approx- 
imation using the EKSO show a correlated structure, 
while differences from the TDSE momentum pair den- 
sities (Fig. [I]) remain. 

Using the PP approximation we obtain momentum 
distributions which yield fundamental insight into the 
double ionization processes. However, this still requires 
knowledge of the exact g xc (xi, X2,t) at time t = T af- 
ter the laser pulse, i.e., of the exact pair density in 
real space. Approximating g xc (x\ ,X2,t) is a formidable 
task itself. This can be seen from the highly correlated 
structure in Fig. [5] where contour plots of the exchange- 
correlation function g yLC {x\ 1 X2, T) are shown for intensi- 
ties where NSDI dominates. An adiabatic approximation 
using the groundstate pair density [8| is not feasible as 
the exchange-correlation function in the entire „4(He 2+ ) 
is required in Eq. (I10| . An expansion for small inter- 
electron distances [22l [23| will not include the correla- 
tions for large \x±— a^l, which are clearly present in Fig. [5] 
By multiplying the complex exchange-correlation func- 
tion with a damping function F(\xi — a^l) with F — ► 
for large |xi — X2I, we verified that short-range correla- 
tions alone in the final wavefunction do not reveal the 
characteristic peaks in the ion momentum density. It is 
therefore of central importance to devise new strategies 
of approximating g xc (xi , X2, t) . 



VI. SUMMARY 

A model Helium atom in strong linearly polarized few- 
cycle laser pulses was investigated. Solution of the time- 
dependent Schrodinger equation yielded momentum pair 
distributions of the electrons freed in double ionization 
and corresponding ion momentum densities. They were 
consistent with a recollision process and, at higher laser 
intensities, with sequential double ionization. These re- 
sults served as a reference for a time-dependent density- 
functional treatment of the system. It was shown that 
the choice of the correlation potential in the Kohn- 
Sham equations is of minor importance compared to the 
form of the functionals for calculating the momentum 
distributions. An uncorrelated approach was found to 
produce ion momentum densities differing significantly 
from the Schrodinger solution in qualitative terms. We 
constructed an exact correlated functional via the two- 
electron wavefunction. The product-phase approxima- 
tion reduces the problem of approximating this func- 
tional. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. 
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